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We investigate the existence of brane solutions in braneworld scenarios described by real scalar 
field in the presence of modified f(R) gravity with a single extra dimension. In the case of constant 
curvature, we obtain first-order differential equations which solve the equations of motion and ease 
the search for explicit analytical solutions. Several examples of current interest are investigated to 
illustrate the results of the present work. 
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1. Introduction. The brane theory is being intensively investigated as a sound candidate for solving the 
hierarchy problem, and other fundamental problems in high energy physics [1-5] . 

The brane concept which was first introduced in [1] is formulated in terms of a single infinite extra dimension 
according to which the physical world appears as a four-dimensional spacetime embedded into an anti-de Sitter 
(AdS) space of one higher spatial dimension. This scenario starts with the five dimensional spacetime AdS*,, 
with warped geometry described by a single function which is real and depends only on the fifth or bulk 
coordinate, the extra dimension, which is an infinite spatial coordinate. The scenario is admissible under the 
presence of scalar fields which are also supposed to depend only on the extra dimension and work to stabilize 
the braneworld configuration [1-6]. We name this the standard braneworld scenario, since it is governed by 
some Einstein-Hilbert action which includes gravity and scalar fields in a standard way. Similar ideas can also 
■ be found in [7]. 

The development of the brane concept was done almost simultaneously with other investigations, initiated 
t-H ■ with the observations of [8] , which has led to the modern concept of dark energy, a generic name used to identify 
the need of an extra but very important component of energy density to fill the spacetime and give rise to the 
present acceleration of the Universe - detailed investigations can be found in the recent revisions [9] . 

We recognize as the standard cosmological model the ACDM model, which describes standard gravity with 
the scalar curvature R, with the inclusion of the cosmological constant A and the nonrelativistic dustlike matter. 
The presence of dark energy may be considered with diverse modifications of this standard scenario. The two 
f---^ main directions are the quintessence way, in which one considers the possibility of changing the standard scenario 
with the inclusion of scalar fields, and the f(R) way, in which one modifies gravity, changing R — ► f(R), where 
J> . f(R) stands for a given function of the scalar curvature. 

The recent interest on f(R) gravity models [10] has increased very radically, and this has inspired us consider 
the issue of investigating the f(R) modification of the braneworld scenario implemented in [1]. The subject of the 
present work is then a natural idea, which consists in the search of brane solutions with the f(R) modification 
of gravity. Even knowing that this problem is hard to solve on general grounds, we believe that our recent 
advances on the first-order framework [11] applied to study scalar fields coupled to gravity may help us to find 
explicit solutions. 

Other ideas of modifying the scenario of [1] have appeared in [12, 13], with the change of the scalar field 
dynamics to the tachyonic dynamics [12] or the inclusion of a Gauss-Bonnet term [13]. Although our proposal 
of including a f(R) modification of gravity is different, it will certainly contribute to the subject. 

2. The framework. The standard braneworld scenario is described by real scalar field <f> interacting with 
gravity via the usual Einstein-Hilbert action, which has the general form 

d i xdy^(--R + c] (1) 
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where we are using standard notation, and taking AttG — 1 . The scalar field <j) is governed by standard dynamics 
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with Lagrange density given by 



C= l -g ab d a <pd b <i ) -V^) (2) 



where g a b describes the five dimensional spacetime, with a, b = 0, 1,2,3,4 and x 4 = y standing for the extra 
dimension. The line element is given by 

ds 2 = e 2A rj^dx^dx" - dy 2 (3) 

where e 2A is the warp factor, and A — A(y) a real function of the extra dimension which gives rise to the warped 

geometry and rj^ = diag{-\ ) describes the four-dimensional flat spacetime, with fj,, v = 0,1,2,3. The 

geometry of the five dimensional spacetime is then described by A(y), and is driven by the extra coordinate y 
alone. We name this the standard flat brane scenario. 

In this work we maintain the metric (3) and scalar field dynamics (2), but we consider the f(R) modification 
of gravity. This means that in the above action we change R — > f(R) to get to the modified action 

S' = f d 4 xdy^g (-\f{R) + Cj (4) 

which will be the action to govern the modified braneworld scenarios to be considered below. 

The five dimensional metric given above implies that the scalar of curvature depends only on the fifth coor- 
dinate y, that is, R = R(y). It can be written in terms of A = A(y) in the form 

R(y) = 8A" + 20A' 2 , (5) 

where the prime denotes derivative with respect to y. 

We consider the modified action (4) and suppose that the scalar field only depends on the extra dimension, 
that is, = 4>{y). Variation of the modified action with respect to the scalar field leads to the equation of motion 

^ + 4AV-^=0 (6) 

dep 

which is the standard equation. However, the f(R) modification changes the Einstein equations and we now 
get to 

RabfR - \gabf{R) + (9abV 2 - V a V 6 )/fl - 2 T ab , (7) 

where fa represents the R-derivative of /, that is, fn = df(R)/dR. We are also using 



2 Sjy/gC) 
V9 Sg a 



T «b = — c ab (») 



as the energy- momentum tensor corresponding to the C given by Eq. (2). 

The Christoffcl symbols associated to the metric (3), have non zero components = A 1 do and T% — 
r}f} 1 e 2A A' . As the A field depends only on y, the components of the Ricci tensor are R^ v — r\ ilv e 2A {A" + 4A' 2 ), 
R44 = —4(A" + A' 2 ), and R^ = R^a = 0. Thus, the Ricci scalar takes the value given in Eq. (5). Also, the 
components of the energy-momentum tensor are given by T^ v = r] lll ,e 2A (^(j)' 2 + V((f))), T 44 = j<f>' 2 — V{<j>) and 
= T^4 = 0. 

We substitute the above results into (7) to obtain the two equations 

\(A" + AA' 2 )f R l -f{R) {\A'f' R + = \<t>' 2 + \V{4>) (9a) 

{A" + A' 2 )f R \f{R) A'f' R = -\cp 12 + iv(4>) (9b) 
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or better 

A"f R -{\A'f' R -\&) = (10a) 
{A" + A' 2 )f R ±f(R) A'f R = + \v{<j>) (10b) 

These are the Einstein's equations, written appropriately to the present study; if we take f(R) — R and use (5) 
they become 

A" = -\^ and A' 2 = -\V{4>) (11) 

which are just the equations we usually investigate in the standard scenario [5]. 

The above equations (6) and (10) are the equations we have to deal with to find solutions in the modified 
braneworld scenario. In the standard scenario, we know that there are first-order differential equations which can 
be applied to solve the equations (6) and (11). Thus, in the next section we consider some specific possibilities 
of writing first-order equations in the modified case, governed by (6) and (10). 

3. The constant curvature case. Let us first consider the constant curvature case. We have that R' = 
and this reduces the equations (10) to the form 

3f R A" + 2(f)' 2 = (12a) 
2A' 2 f R \f{R) = V(J>) (12b) 

In first place we describe the case of an empty space. We consider <j> = <j> = const, with V{4>) = — 1/L 2 , for L 
being a lenght scale. These equations are thus reduced to 

3f R A" = (13a) 
2A' 2 f R - - A f{R) + -L = (13b) 

If fa 7^ 0, we find that A" = 0; thus, the brane solution is A(y) = —\y\/L; the singularity at the brane location 
is the same of the standard case sec, e.g., [3]. Thus, R — 20 /L 2 . As a result we get the following equation for 
the curvature 

k i n 

RfR--f(R) + jz=0. (14) 

This equation can be treated in two ways: (i) we suppose that f(R) = R+"/R n in order to arrive at the equation 
(n — 5/2)jR n — (3/2)R + 10/i 2 =0, from which we can find R; (ii) we can treate this equation as a differential 
equation for f(R), whose solution is 

f(R) = ^ + CR b ' 2 . (15) 

In a more general case, let us now consider R constant in the presence of the scalar field. We thus implement a 
first-order framework to solve the system. It is inspired in the investigations [11], and is given by the suggestion 
that 

A'{y) = W(<j>(y)). (16) 

This establish an explicit relationship between A and <f) via the function W = W{4>). Then, using the equations 
above we can write 

<t>'(y) = -\fn w* (17) 



4 



This pair of first-order equations allows to obtain 



and the potential 



A" = -\f R Wl (18) 



V{4>) = -\f{R) + Vr w 2 ^f R w 2 (19) 



In principle, W is assumed to be an arbitrary function. However, the form of the metric (3) restricts it 
severely. In fact, using the ansatz (16), the scalar curvature R can be written as 

R = 8A" + 20 A' 2 = -12 f R W 2 + 20W 2 (20) 

Therefore, W has to obey the constraint 

w *-i; w2+ ik R=0 (21) 

It shows that the form of W depends explicitly on the specific scenario, that is, on the specific modification 
introduced by f(R). 

If we now use (17) and (21) we can write 

= -'lfmj + c, (22) 



/■ #. = [ 

J ™* J M 



VPiW 2 - A 
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where we have set (3\ — 5/ (3/i{), /?2 = R/ (12/i?) and c is an integration constant. We see that for obtaining the 
field <p(y) we need to first calculate W((f>), and then integrate and invert Eq. (22), which may make the problem 
very hard to solve. 

To study Eq. (21) it is convenient to write it as 

^|-/W 2 +/3 2 = 0, (23) 

which allows us identify the following characteristic cases: 

1. Zero curvature (R = 0). 

(a) Case f R > (or (3\ > and fii =0). In this case we have to solve 

This equation can be integrated to give 

W(<f>)=exp(± y /^(<f,-b 1 )), (25) 
and the corresponding field solution reads 

4>{y)=T-LMT{5/2y-cfh)) + c 1 . (26) 
VPi 

Here b\ and c\ are integration constants. 

(b) Case Jr<0 (or /3i < and f3 2 = 0). This case presents no real solutions, as f3\ < in equation (24) 
appears inside a square root. 

To obtain the field A(y) we can directly solve the equation R = 8A" + 20A' 2 = 0, which yields 

e 2A(y) = ai {y-y,V (27) 

with a\ and j/q being integration constants. 



2. Negative curvature (R < 0). 

(a) f R > (i.e. /?i > and /3 2 < 0). We have 



— = ±V/W 2 + I/? 2 |. (28) 
Integration of this equation yields 

W(0) =±2/3sinh(yA"0 + 6 2 ), (29) 

and the field solution for this case is 

<P{y) = ±^(arcsmh(tan(/3(5y-2c/3i))) + c 2 ). (30) 
vPi 

Here & 2 and c 2 are integration constants, and (3 = 5\/|/3 2 |/|/3i|- 

(b) The case f R < (or /?i < and /3 2 > 0), is peculiar. Since W(<p) turns out to be purely imaginary 
this case is not allowed. 

For the negative curvature case we can write the equation for A(y) as R = 8A" + 20A' 2 = — \R\- The 
solution yields 

e2Mv) = (f£f (™(lVm(y-yo))y (3i) 

where a 2 is an integration constant. 
3. Positive curvature (R > 0). 

(a) f R > (i.e. p 1 > and /3 2 > 0). We have 



(h> ±V/W 2 -/?2- (32) 
This equation is easily integrated to give 

W{4>) = ±2/?cosh(V/3i> + &3) ; (33) 

with & 3 being an integration constant. In this case, integrating and inverting equation (22) leads to 
inconsistent solutions for <p{y). 

(b) f R < (or /?i < and (3 2 < 0). We have 



— = ±Vl/?2|-|/3i|^ 2 , (34) 
which gives the solution 

W(<£) = ±2/?sin( v ^]0 + &4), (35) 
and the corresponding field solution 

4>{y) =±^^(arcsin(tanh(/3(5y + 2c|/3 1 |)))+c 4 ), (36) 
VlPil 

with 64 and C4 being integration constants. 
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When R is positive, we can find the warp factor by solving R = 8 A" + 20 A' 2 . The solution for this 
equation looks like 



e 2A(y) 



(f) 5 ( cosh Q^(y-yo))) 5 , (37) 



with ci4 being an integration constant. 

In the standard scenario, the presence of first-order equations is directly connected with supersymmetry and 
they nicely appear within the supergravity context [3, 7]. In this sense, the first-order framework which we have 
just obtained suggests that the model engenders supersymmetric extension. 

4. Other procedures. The first-order framework is very interesting, since it helps in the search of analytical 
solutions. However, its application in the form used in the previous section does not exhaust all possibilities. We 
then make further efforts below. In particular, we may propose specific ansatzes as a direct way to investigate 
the problem. In this case, it is better to consider the system of equations given by Eqs. (6) and (9a), which we 
rename as 

^ + 4A'<j>> ^1 = (38a) 
2(A" + AA' 2 )f R - f(R) - ( 6A'f' R + 2/£) - 20' 2 + 4U(0) (38b) 

We study this system for the specific choices which follow below. 

4.1 The case of constant curvature. If R is constant, the equation (38b) is reduced to 

-/OR) + 2 (4A' 2 + A") f R (R) = 20' 2 + 4U(0). (39) 

As the function A(y) is determined by the condition on the scalar curvature R, it remains to find <fi(y), V(<j>) 
and f(R) for all the three possibilities, R < 0, R — 0, and R > 0. 

Let us first consider the case A(y) = — k\y\. For this A, we have R = 20k 2 and A" = for y ^ 0. Thus, after 
differentiating eq. (39) we can eliminate the potential term using (38a) to obtain an ODE for <p 

(/>'(<(>" + 2A'<f>') =0. 

Using that A' = k sgn(y) we obtain the solution 

</> = 0o e~ 2%l + Oo , (40) 

with 0o and constants. Therefore, from equation (38a), the consistent potential is 

V(4) --2fc 2 (0-Oo) 2 . (41) 

Because of the form of A, the right-hand side of (39) turns out to be zero, providing an equation for f(R). 
Using that A(y) = —k\y\ and R — 20k 2 , this equation takes the form 

f(R) - ^Rf R = 0, (42) 

which has the solution 

f(R) = aR 5 / 2 . (43) 



with a = const. Thus, we have found another solution for the constant positive curvature case. 
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We can also evaluate the zero curvature case. We study the case f(R) = R + jR n , so f(R) — 0, but fa = 1 
for n > 1. The function A(y) is given by (27) but, for simplicity, we choose y a = using the translational 
invariance. Therefore, we must solve the system 



1_ , _ dg(fl 



^" + 4^-^ = 0; (44) 



2 ^ 2 + 4 ^) = -(|^- ( 45 ) 

A natural ansatz in this case is 

cf>(y) = kln\y\, V (<]>) = b e a * . (46) 

Substituting this ansatz we find that the values of k, a and 6 are fixed as k = 3±-\/3, a = —2 fc _1 , and 6 = — 3fc/5, 
while 7 remains a free parameter. 

There exists another solution for this zero curvature case (in which, however, it is difficult to express the 
potential Vasa function of <f>) 

<t>(y) = ±VTE (-l + ayi)?+ arctan((-l + ciy^)"^ + c 2 , (47) 

Here ci and c 2 are integration constants. 

4.2. The case of non constant curvature. We can also try to solve the system (38) without restricting 
the curvature to a constant. For the function f(R) in the form f(R) = R + 7-R™ the system becomes 

^ + 4A>'-^M = 0; (49a) 
a<p 

-R - ~/R n - 2(j)' 2 - W{4>) + [8A' 2 + 2A")(1 + n 7 i?"- 1 ) - 

-6n(n - l)jA'R n - 2 R' + 2 [n(n - l)(n - 2)jR n - 3 R' 2 + n(n - l)R n ~ 2 R"] = 0. (49b) 
Here we can substitute the curvature (5) together with the expressions 

J?' = 84'" + 4CL4/A", 

R" = 8 A"" + AO A" 2 + AO A' A 1 " . (50) 

Then the system (49) reads 

^ + 4Ay-^M=0; (51a) 

-8A" - 20A' 2 - 7 (8A" + 20A' 2 ) n -2<f>' 2 - W(<j>) + 
+(8A' 2 + 2A")(1 + n- f {8A" + 20A 1 2 ) n - 1 ) - 
-6A'n(n - l)j(8A" + 20A' 2 ) n - 2 [8A"' + AO A' A"] - 
- 2n(n - 1) 7 [(n - 2)(8A" + 20A' 2 ) n - 3 (8A'" + A0A l A") 2 + 

+ (8A" + 20A' 2 ) n - 2 {8A"" + AOA" 2 + AO A' A 1 ")} = 0. (51b) 

This system can be solved using a variation of the first-order formalism used before. Let us suggest the 
following choices of ansatzes 

A' = Ce™*, <t>' = De m *. (52) 
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This is convenient since it allows to decouple the first equation (51a) which is now treated as a differential 
equation for the potential. We note that it follows naturally from the above choices that A = A</>, with A = C/D 
being a constant. One can find A" = CDme 2 "^, A"' = 2CD 2 m 2 e 3m <t> , A"" = 6CD 3 m 3 e im ' t >, <j>" = D 2 me 2m f 
Substituting these expressions into (51a) leads to 



which is easily integrated to give 



dV 

(D 2 m + ACD)e 2m * - — = 0, (53) 
d<f> 



VU) = — (mD 2 + 4CT>)e 2m< ^, (54) 
2m 



which has an exponential behavior, of current interest to cosmology - see e.g., [9]. 

We substitute the above choices (52) into the second Eq. (51b), and we use the exponential potential (54) to 
get to the following equation 



e 2m<p 



(SCDm - 20C 2 ) - 7 (8CT>m + 20C 2 )"e 2nm * - (ad 2 + ^-^j e 2m<i> + 

+ (8C* 2 + 2CHm)[e 2m * + n 7 (8CHm + 2OC* 2 )"" 1 e 2nm0 ] - 

-6Ce 2nm ^(n - 1) 7 (8CT>to + 20C 2 )™- 2 (16CT> 2 m 2 + 40C 2 Dm) - 

-2n(n - l) 7 e 2nm< *[(n - 2)(8CDm + 20C ,2 )"- 3 (16CT> 2 m 2 + 40C 2 Dm) 2 + 

+ (8CDm + 20C 2 )"- 2 (48CT> 3 m 3 + 120C 2 L> 2 to 2 )] = 0. (55) 

Comparing the terms associated with e 2m * and e 2nrn ^ (recall that 1) and extracting the overall factors D 2 
and D 2n , we find two algebraic equations for to and X — C/D 

(m + 2A)(3Am + 2) =0; 

-(8Am + 20A 2 )™ + n(8A 2 + 2Am)(8Am + 20A 2 )™" 1 - 
-6\n(n - l)(8Am + 20A 2 )"- 2 (16Am 2 + 40A 2 m) - 
-2n(n - 1) [{n - 2)(8Am + 20A 2 )"~ 3 (16Ato 2 + 40A 2 m) 2 + 

+ (8Ato + 20A 2 )"~ 2 (48Ato 3 + 120A 2 m 2 )] = 0. (56) 

We see that there are two possible solutions for A in terms of m: A = —2/ (3m) and — to/2. Below we find a way 
to fix A as a function of to. The corresponding solutions for the fields (f> and A following the choices (52) are 
given by 

4>{y) = -—ln{C -mDy), 

TO 

A(y) = --HCo-mDy), (57) 

TO 

where Co is an integration constant. The curvature in this case is 

The singularity shown at j/o = Cq/uiD identifies the location of the brane. Of course, the singularities that we 
have found in the present work are not specific of the f(R) modification here studied - similar singularities also 
appear in the standard scenario - see, e.g., Rcf. [14]. 
We now use f(R) = R + jR n in Eq. (10a) to get 

A" {I + n-fR 11 - 1 ) - \ in {n - l)R n ~ 3 [A'R'R - (n - 2){R') 2 - R"] = -\(j>' 2 . (59) 
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This equation and the solutions (57) lead to another equation, which can be written with two kind of terms, 
one proportional to 1/(C — mDy) 2 and the other to 1/(C — mDy) 2n : 

Ti T 2 
(Co - mDy) 2 + (Co - mDy) 2n = °" ^ 

Thus, since n ^ 1, from T\ we obtain A = —2/ (3m), and so the second possibility A = — m/2 has to be discarded. 
We then use A = -2/ (3m) in T 2 = to get 

"■^^(rb-i)- (81) 

which leads to the possible values of m, for several 1. 

5. Ending comments. In this work we have extended the standard braneworld scenario to include the f(R) 
modification of gravity. This modification makes the problem very hard to solve, but in the case of constant 
curvature we have been able to find a pair of first-order equations which solve the second-order equations of 
motion, and we have investigated several models. The presence of first-order equations eases the investigation, 
and we have included explicit solutions for all the three cases of positive, null or negative curvature. 

As we know, in the standard scenario the presence of first-order equations is directly related to supersymmetry. 
Thus, the first-order framework here obtained suggests that the f(R) modification introduced in the present 
work may engender supersymmetric extension. Another issue concerns stability of the investigated scenarios, 
with the f(R) modification here introduced, which may be studied within the lines of the third work in [11]. 
These and other issues are presently under consideration, and we hope to report on them in the near future. 
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